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I.        INTRODUCTION 
A.       BACKGROUND 

The  use  of  a  natural  medium  for  communications  systems,  especially  for  wireless 
and  mobile  radio  communications  systems,  implies  unavoidable  degradation  associated 
with  the  randomness  that  accompanies  natural  phenomena.  Refraction,  reflection, 
diffraction,  scattering,  focusing  attenuation  and  other  factors  cause  variations  between 
the  received  signal  level  and  a  calculated  free-space  level  for  a  given  transmitter  output. 
The  degradation  is  usually  due  to  natural  causes  but  can  also  be  caused  by  man  made 
items  such  as  huge  buildings  in  an  urban  area  [1,2]. 

Although  only  the  direct  wave  is  desired,  several  distinguishable  paths  between  the 
transmitter  and  receiver  do  exist.  These  are  typically  due  to  reflections.  These  paths  cause 
several  waves  to  arrive  at  the  receiver  at  slightly  different  times  and  produce  fading.  The 
interference  between  a  direct  wave  and  a  reflected  wave  is  called  multi-path  fading  [2]. 

The  electro-magnetic  waves  tend  to  arrive  at  the  receiver  over  different  paths 
having  different  propagation  delays.  The  received  signal  consists  of  components  having 
randomly  distributed  amplitudes,  phases  and  angles  of  arrival.  When  there  is  no  single 
line-of-sight  path  between  the  transmitter  and  the  receiver,  the  received  signal  is  modeled 
as  a  Rayleigh  random  variable.  Its  probability  density  function  (pdf)  is  given  by  [4] 


f(z)= 


z 

—  exp 


(      z    ^ 


2(7  2 


0  otherwise 


z>0 

(1) 


where  2<72is  the  average  power  of  faded  signal,  and  z  is  the  amplitude  of  the  received 
signal. 


The  resulting  communication  channel  is  called  a  Rayleigh  fading  channel.  When  the  line- 
of-sight  component  is  not  zero,  the  received  signal  is  modeled  as  a  Ricean  random 
variable.  Its  pdf  is  given  by  [4] 
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where  a2 is  the  average  power  of  the  unfaded  (direct)  signal  component,   2<72is  the 
average  power  of  the  faded  signal  component  and  70(«)  is  the  zeroth-order  modified 

Bessel  function  of  the  first  kind.  This  resulting  communication  channel  is  called  a  Ricean 
fading  channel  [2,3]. 

Rayleigh  fading  severely  degrades  the  average  probability  of  the  bit  error 
performance  of  the  receiver.  One  way  of  minimizing  this  problem  is  to  increase  the 
transmitter  power  or  the  antenna  size.  However,  such  a  solution  requires  costly 
engineering  alterations  making  other  alternatives  very  attractive.  It  is  well  known,  that 
diversity  combining  is  one  of  the  most  effective  alternatives  [1]. 

To  reduce  the  effects  of  fading,  one  of  several  diversity  techniques  can  be 
implemented.  Diversity  is  a  procedure  to  receive  or  transmit  the  same  symbol  multiple 
times  in  order  to  provide  redundancy  at  the  receiver.  The  basic  idea  of  diversity  is  that 
some  of  the  received  redundant  symbols  will  be  more  reliable  than  the  others,  and  the 
demodulation  decision  will  be  made  using  the  more  reliable  symbols.  In  order  to  be 
useful,  each  redundant  symbol  must  be  received  independently  [4].  The  diversity  can  be 
implemented  as  space,  time,  frequency,  angle  and  polarization  diversity.  Space  diversity 
consists  of  using  multiple  antennas  at  the  receiver  in  order  to  receive  the  transmitted 
symbol  multiple  times  via  multiple  paths.  Time  diversity  is  employed  by  transmitting  the 


same  symbol  multiple  times.  In  frequency  diversity,  the  symbol  is  transmitted  on  multiple 
carrier  frequencies  at  the  same  time.  Angle  diversity  is  implemented  by  a  set  of  directive 
antennas  where  each  one  responds  independently  to  a  wave  to  produce  an  uncorrelated 
faded  signal.  Polarization  diversity  is  possible  with  two  orthogonal  polarizations. 
Combining  the  signals  from  the  diversity  branches  allows  a  significant  performance 
improvement.  Diversity  combining  can  be  accomplished  prior  to  or  after  detection,  which 
are  called  pre-  and  post-detection,  respectively.  Diversity  combining  can  be  linear  or  non- 
linear. In  this  thesis,  only  linear  combining  methods  are  analyzed.  Unlike  the  other 
methods  used  for  reducing  fading  effects,  diversity  combining  lowers  costs  and  reduces 
complexity  in  the  receiver  [2].  In  this  thesis,  a  new  method  referred  to  as  Post-Detection 
Selection  Combining  (PDSC)  is  presented  and  compared  with  Selection  Combining  and 
Equal  Gain  Combining  [1,3]. 

Post-Detection  Selection  Combining  (PDSC)  is  invented  by  Professors  Tri  T. 
Ha  and  Ralph  D.  Hippenstiel  [14].  This  method  selects  the  signal  with  the  largest 
amplitude  or  combines  the  two  (three)  received  signals  with  the  two  (three)  largest 
amplitudes  after  detection.  Since  combining  is  done  after  detection,  the  new  method  is 
named  as  First,  Second,  and  Third  Order  Post-Detection  Selection  Combining  and  is 
denoted  by  PDSC-1,  PDSC-2  and  PDSC-3,  respectively. 

B.        SYSTEM  DESCRIPTION 

A  noncoherent  differential  phase  shift  keyed  (DPSK)  receiver,  as  shown  in  Fig.  1, 
is  used  to  demonstrate  the  performance  of  Post-Detection  Selection  Combining.  As  a 
noncoherent  receiver,  the  DPSK  receiver  requires  no  phase  recovery  and  can  be  build 
inexpensively.  For  this  reason,  it  is  widely  used  in  wireless  communications  systems. 


While  DPSK  signaling  has  the  advantage  of  reduced  receiver  complexity,  its  energy 
efficiency  is  inferior  to  that  of  coherent  PSK  by  about  3  dB  [4].  The  average  probability 
of  error  for  DPSK  in  additive  white  Gaussian  noise  is  given  by  [4] 


Pe=2e: 


(  _ 

(3) 


~Eb 


No 


where  Eb  I N0  is  the  bit-energy-to-noise-density  ratio. 

Before  we  analyze  the  Post-Detection  Selection  Combining,  we  need  to  consider  general 

structure  of  the  DPSK  receiver.  The  phase  difference  between  two  successive  bits  is  used 

to  convey  information.  Whenever  the  data  bit  bt  =  0  is  transmitted,  the  waveform  for  two 

consecutive  differentially  encoded  bits  ci_xci  is  given  by  [6] 

vm(t)  =  ±ApT(t-(i-\)T)cos(2nfct+d)±  A pT  (t-iT) cos (2 n  fc  t  +  d) 

=  [±ApT(t-(i-\)T)±A  pT  (t -iT)]cos(2n  fc  t  +6).  (4) 

When  the  data  bit  bt  =  1  is  transmitted,  the  waveform,  which  has  a  phase  change  of 
k  radians  between  two  consecutive  bits,  is  given  by 

v(2)  (0  =  ±ApT(t-(i-  l)T)  cos  (27C  fct  +  0)+  ApT(t-  iT)  cos  {in  fc  t  +  d) 

=  [±  ApT(t-(i-l)T)+A  pT{t-iT)]cos(2Kfct+9),  (5) 

where  6  is  the  signal  phase,  fc  is  the  carrier  frequency,  A  is  the  amplitude  and  the 

±  signs  denotes  the  polarity  of  the  two  consecutive  bits.  The   pulse  pT  (t  -  iT)  is  given 
by 

pT(t-iT)=[l   '  'T<f*(*+1>rJl  (6) 

[0  ,  otherwise 

When  bx  =  0  is  transmitted,  the  received  signal  is  of  the  form 


4 


r(t)  =  [±ApT(t)±ApT(t-T}]cos(2Kfct  +  e)+n(t),  (7) 

where  n(t)  is  additive  white  Gaussian  noise  with  zero  mean  and  power  spectral  density 


The  in-phase  outputs  of  the  receiver  are  given  by 
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The  quadrature  outputs  of  the  receiver  are  obtained  as 
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(12) 


(13) 


Nlc ,  N2c ,  Nu  and  N2s  are  independent,  identically  distributed  (i.i.d.)  Gaussian  random 


variables  with  zero  mean  and  variance  equal  to 


NJ 


Assuming  that  bl  =  0  is  transmitted,  the  outputs  of  the  signal  branch  are 

Ylc=±AT  cos  6  +  Nlc, 
Yu=+ AT  sind  +  NU, 

and  the  outputs  of  the  non-signal  branch  are 

Y2c  =N2c, 
Y2s  =N2s. 

As  seen  in  Fig.  1 ,  the  decision  variable  for  the  signal  branch  is  given  by 

Yu2+Yu2=(±ATcose  +  Nj2+{+ATsmd  +  Nj2. 

For  the  non-signal  branch  ,  the  corresponding  variable  is  given  by 

Y  2  +Y  2  =  N    2  +N    2 


(14) 
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Figure  1 .  Block  diagram  of  a  noncoherent  DPSK  receiver 


For  convenience  ,  we  normalize  (16)  as  follows: 
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By  setting  F,  =  j(Y2  +Y;s),  nlc  =  JjNlc  .  «,,  =  nfc#i,  .  A  =  ±-^r  A'  we  have  the 


decision  variable 


7,  =  (2A'cosO  +  nu  f  +{-2A'sm6+nu  f 


(19) 


For  a  Rayleigh  fading  channel,  A'  is  a  Rayleigh  random  variable  with  parameter  <T  as 
given   in   (1).   Furthermore,   for  a  uniform  variable    0,   A'cos#   and   -A'sin#    are 

Gaussian   random   variables   with   zero  mean   and   variance    <72[11].   The   variables 
2A'cos0   and  -2A'sin0   are  both  zero  mean  Gaussian  random  variables  with  variance 

4c 2  .The  zero  mean  Gaussian  random  variable  nu  and  nu  are  independent  with  variance 

(J2  =  N0  .Therefore  the  random  variable  Yi  (19)  is  a  central  chi-square  random  variable 


with  the  following  density  function  [4] 
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exp 


2cr,2 
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(20) 


where 


(J2  =Ao2  +<J2  =4<72  +  N0 


(21) 


Define  the  average  energy  per  diversity  channel  as 

E  =  2(J2; 


(22) 


then,  for  L  diversity  channels,  the  bit  energy. is 

Tb=LE.  (23) 

For  the  non-signal  branch,  the  decision  variable  Y2  =  —  \Y2c  +  Y2s )  can  be  written  as 

Y2=n22c+n22s,  (24) 

[2  [2 

where  rt2c  =  J—N2c  and  n2s  =  J—N2s  are  zero  mean  Gaussian  random  variables  with 

variance  a  I  =  N0 .  The  density  function  of  Y2  is  a  central  chi-square  density  function  [4] 

v2>0,  (25) 


A(^)~exp 
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where 


e22=<J2n=N0.  (26) 


C.       OBJECTIVE 

In  the  following  chapters,  probability  density  functions  are  obtained  and  the  bit 
error  probabilities  are  derived  for  PDSC-1,  PDSC-2,  and  PDSC-3.  In  Chapter  VI,  Post- 
Detection  Selection  Combining,  Selection  Combining  and  Equal  Gain  Combining  are 
compared. 

The  primary  objectives  are  to 

(a)  present  the  new  Post-Detection  Selection  Combining  method, 

(b)  make  the  necessary  bit  error  rate  derivations,  and 

(c)  compare  the  performances  of  Post-Detection  Selection  Combining,  Selection 
Combining  and  Equal  Gain  Combining. 


II.       REVIEW  OF  PREVIOUS  WORK 


A.         EQUAL  GAIN  COMBINING 

The  Equal  Gain  Combining  (EGC)  method  adds  in  a  noncoherent  fashion  equally 
weighted  branch  signals.  Although  EGC  is  one  of  the  most  commonly  used  diversity 
techniques,  the  receiver  is  dependent  on  the  order  of  the  diversity  [3].  The  bit  error 
probability  expression  for  Equal  Gain  Combining  is  given  by  [4] 
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(27) 


B.        FIRST  ORDER  SELECTION  COMBINING  (SC-1) 

In  first  order  selection  combining  (SC-1)  the  signal  with  the  largest  amplitude 
(hence  the  largest  signal-to-noise  ratio)  is  selected  as  shown  Fig.  2. 
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Figure  2.  Selection  Combining  with  L-th  order  Diversity 


The  conditional  bit  error  probability  for  SC-1  is  simply  the  bit  error  probability  for  DPSK 
conditional  on  the  signal-to-noise  ratio  y  [4] 


pb(r)=\e-y 


The  pdf  of  y  is  given  by  [3] 
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where 
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(30) 


Here  E  is  the  average  energy  per  diversity  bit,  and  El N0  is  the  average  diversity  bit 

energy-to-noise  density  ratio. 

The  bit  error  probability  for  SC-1  can  be  evaluated  as  follows  [3] 


Pb=]pb(r)f(r)dr- 


(31) 


Inserting  (28)  and  (29)  into  (31)  and  performing  the  integration,  we  obtain 
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We  note  that  the  average  bit  energy  Eb  is  related  to  the  average  diversity  bit  energy 

_        _  ln  ~Eb 

E  by  E.  =  LE .  Thus  using  the  identity  a  =  _  °  ,  we  obtain  a  result  in  terms  of as 

Eb  N0 


fLN^ 


L-\ 


L  k=0 


<    k    > 


(-If 


1  + 


'LNn^ 


k  Eb 


+  k 


K  Eb    , 


(33) 
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C.        SECOND  ORDER  SELECTION  COMBINING  (SC-2) 


In  second  order  selection  combining  (SC-2),  two  signals  with  the  largest 
amplitudes  or  signal-to-noise  ratios  at  the  input  of  the  DPSK  demodulator  are  combined. 
The  bit  error  probability  conditional  on  the  signal-to-noise  ratio  y  is  given  by  [3] 


A(r)=^-r[4+y]- 

o 


The  pdf  of  y  is  also  given  in  [3] 


f(y)=L{L-l)ae-^ 


L-2 


(L-2)(-l)k 
k 


f  _aky  \ 

l-e~  2 


\  ) 


/J 


(34) 


(35) 


where 


E 


as  defined  in  (30).  Substituting  (34)  and  (35)  into  (31)  yields 


(36) 


Pb=j^[4  +  y]L(L-\)ae-ar 


2  +^ 


*=i 


ayk  \ 

\-e    2 

JA 


dy  .  (37) 


By  performing  this  integration,  the  bit  error  probability  is  found  to  be  [3] 
(  a    Y  L(I-1) 


Pb  = 


v1+«y 


8 


L-2 


xl2  + 


— +2 


L-2 


(-If 


(l  +  a)[4(2  +  2cc  +  ak)+4]+ak 
(2  +  2a+akf 


(38) 
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The  bit  error  probability  is  obtained  in  terms  of as 


Pb  = 


(    LN0     } 


LNC 


L(L-l) 


2  +  - 


LNQ 


L-2 


+ 


(L-2\ 


jt=i 


k 


{-if 


f      LNA  (         LNn     LNn  A 


1  + 


A. 


2  +  2^±  +  ^±k 
Eb        Eb 


+  4 


LNa 
+  =±k 


(         LN0     LNa     V 
K  Eb        Eb 


(39) 


D.         THIRD  ORDER  SELECTION  COMBINING  (SC-3) 

For  the  third  order  selection  combining  (SC-3),  three  signals  having  the  three 
largest  amplitudes  or  signal-to-noise  ratios  at  the  input  of  the  DPSK  receiver  are 
combined  and  the  conditional  bit  error  probability  for  DPSK  is  given  by  [3] 


Pb(r)=— e 


32 


I6  +  67  +  -72 

2 


(40) 


The  pdf  of  7  for  SC-3  is  given  by  [3] 


/(y)=L(L-lXL-2)gg_cr 


X- 


a2_f_    ^(L-^-tf 


k=\ 


\  J 


( 


ayk-3 


1-exp 


kay 


W 


JJJ 


(41) 


Substituting  (40)  and  (41)  in  (31)  yields 
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b         J   09 


32 


x 


I6  +  67  +  -72 

2 


a2r2     «|rL_3N(_1)* 


Ul-\\l-i) 


ae 


-ay 


■+2 


* 


it2 


&a/-3 


1-e"  3 


dy. 


(42) 


The  bit  error  probability  becomes  [3] 


L{L-\\L-2) 


64 


16 


LNf 


fEb^2 


LNr 


LNf 


+  1 


V     '* 


■  +  ■ 


V    f 


—       \5 


LNC 


+  1 


1  + 


7.-3 


1 6k 


48 


l+Eb 


\2k 


'10 

LN„ 


LV 


18 


LNr 


J 


\+T 


LNr 


144 


1  + 


V 


LNr 


3  +  k  +  3 


\ 


LNr 


LNr 


162 


+• 


1  + 


Eb 


■+■ 


LNr 


3k 


Lb 


1+ 


E^2     < 


LNr, 


3  +  fc  +  3 


Eb 


LN 


11 


LNr 


1  + 


LNr 


LNr 


81 


<E>* 


LNr 


1+ 


r^ 


LNr 


3  +  k  +  3 


V 


LNr 


<43) 


In  the  next  chapter  we  derive  the  bit  error  probability  when  first  order  post- 
detection  selection  combining  is  employed. 
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III.     FIRST  ORDER  POST-DETECTION  SELECTION 
COMBINING  (PDSC-1)  ANALYSIS 


In  the  first-order  Post-Detection  Selection  Combining  (PDSC-1)  technique  [14], 
the  maximum  of  L  diversity  samples  at  the  output  of  the  L  noncoherent  DPSK 
demodulators  are  considered.  The  diversity  receiver  is  shown  in  Fig.3. 


r(t) 


cos(27tfct) 


t  =  T,2T,. 


►<gh> 


/(•)*  -^r^  T 


sin(2#"f) 


+ 


+ 


©— ►    (•) 


I 
+f 


max{Yu,Yi,2,....YiL} 


V, 


Figure  3.  Block  diagram  of  a  noncoherent  DPSK  receiver 

with  first  order  Post  -  Detection  Selection  Combining 
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A.       PROBABILITY  DENSITY  FUNCTIONS  OF  DECISION  VARIABLES 

As  shown  in  Fig.  3  V,  and  V2  denote  the  decision  variables  for  the  signal  branch 


and  non-signal  branch,  respectively.  We  have 

V^rnax^,,,^....^}, 

V2  =max{y21,F22....r2X}, 


(44) 
(45) 


where  Ylk  and  Y2k ,  k  =  1,2,... L  are  the  amplitudes  at  the  output  of  the  DPSK  demodulator 

of  the  k-th  diversity  channel.  We  assume  that  Yu  and  Y2k  are  independent  and  identically 

distributed  random  variables  for  all  k. 

For  the  signal  branch,  the  probability  density  function  of  the  k-th   signal 
amplitude  Yu  is  a  central  chi-square  density  function  given  in  (20) 


/",.(**) —exp 


(       y       ^ 

20-' 


7,^0, 


(46) 


where  a2  is  given  in  (2 1 )  and  (22) 


c712=4o-2+A^0=2£:  +  ^0  . 


(47) 


The  probabilty  density  function  of  the  decision  variable    V,    in   (44)  is  given  by 


(Appendix  A) 


/v1(v1)=L/ru(v1)[F,u(v1f-1, 


(48) 


where 


^>W 


o  2cjr 


exp 


v     A 
2cr; 


dyn  =l-exp 


v, 


2c7,2 


(49) 


Substituting  (46)  and  (49)  into  (48)  yields 
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fvM)= 


(   „  \ 


( 


2a; 


exp 


2a; 


1-exp 


/  L 


2c; 


v-li-l 


yj 


(50) 


Similarly,  for  the  non-signal  branch,  the  probability  density  function  of  the  k-th  noise 
amplitude  Y2k  is  a  central  chi-square  density  function  given  in  (25) 


/r2t(3;2it)  =  T-Texp 

2  2(7, 


( 


y 


2  k 


2o; 


y2k^o, 


(51) 


where   o2  =  cr„"  =  Af0  .  The  probability  density  function  of   V2    in  (45)  is  given  by 


(Appendix  A) 


fv2M=LfYJv2)[FY2k(v2)}-\ 


where 


7    1         f    y     }  (     "     ^ 

F*2t  (va)=  j^ — TexP  ~-rk  fyik  =l-exp 


02C72 


2<7 : 


2<xJ, 


(52) 


(53) 


Substituting  (51)  and  (53)  into  (52)  we  have 


/v>2)  =  :r^exP 


f      „     \ 


2<r: 


arj  , 


r    „    \ 


-.L-1 


1-exp 


2cr; 


/J 


(54) 


B.        BIT  ERROR  PROBABILITY 

Using  the  results  derived  in  Part  A,  the  probability  of  error  for  first  order  Post- 
Detection  Selection  Combining  is 


'.-J 


\fvSvi)dvi 


A  (vi )  dvx 


(55) 


Substituting  (50)  and  (54)  into  (55)  yields 
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Substituting  (50)  and  (54)  into  (55)  yields 


'.-J 


2a: 


exp 


2<J22J 


1-exp 


(         ^^ 

v. 


L-\ 


2a: 


dv~ 


x 


2a; 


exp 


V 


2a: 


1-exp 


v^M 


2(7, 


dvx  . 


(56) 


Applying  the  binomial  theorem, 


(*+*)"  =X 


n(jj\ 


Jt=0 


v*y 


(57) 


to  the  exponential  terms  in  (56),  the  exponential  terms  can  be  simplified  as 


-1Z.-1 


1-exp 


2a: 


yj 


L-ifr  _  A 

=X  ,    (-D 


*=0 


V*      J 


(        ( 
exp 


Yf 


2a: 


(58) 


and 


1-exp 
Then,  (56)  becomes 


v, 


lL-1 


2(7,2 


yj 


L-l 


=  1 

7=0 


^L-0 


V    •>        / 


(-iy 


'  r  v  ^ 

exp 


2(7, 


(59) 


*=J 


L    i^fL-O 


ZC7,    i=0 


(-1)' 


f    (     r1+k^ 


V  / 


exp 


-v. 


2a: 


\        "   J 


dv. 


L-l 


2(7,    ;=o 


(L-\\ 


K  J      J 


(-1)' 


r1+Jx<\ 


exp 


-v, 


2(7,' 


V     -'    JJ 


dv, 


(60) 


be 


Performing  the  integration  in  (60),  the  bit  error  probability  for  PDSC-1  is  found  to 


P=tli  L    ^    (L-1)!L("1);" 


k=Q  j=0 


.k  +  lJ 


(L-l-j)\j\(l+j)  i+k  a; 
1+7  <7J 


(61) 


+  1 
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Rewriting  the  parameters  <7,  and  <T2  as 


c}=4<t2+N0=2E  +  Nq 


(62) 


and 


(j22=a2n=N0  , 


(63) 


then,  the  bit  error  probability  is  found  in  terms  of  the  energy  per  diversity  channel  as 


L-\  L-\  f     J       \ 

*=o  j=o  y K  ^ v  j 


(  L   ) 


(-ir 


(     «Vu^ 


1  + 


1  + 


IE 

Nr 


\  +  k 

,  1+7  , 


(64) 


—        Efc 
Substituting  E  by  in  (64),  the  bit  error  probability  for  PDSC-1  is  found  to  be 

Li 


L-\  L-\  (    J       \ 

n=zzLti 


(  L   ^ 
7+1 


(-ir 


i+ 


i+ 


2Eb 


V1  +  ^ 


(65) 


,1+7, 


In  the  next  chapter  we  will  derive  the  bit  error  probability  when  second  order  Post- 
Detection  Selection  Combining  is  employed. 
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IV.      SECOND  ORDER  POST-DETECTION  SELECTION 
COMBINING  (PDSC-2)  ANALYSIS 

In  the  second-order  Post-Detection  Selection  Combining  (PDSC-2)  technique,  the 

two  largest  values  of  L  diversity  samples  at  the  outputs  of  the  L  noncoherent  DPSK 

demodulators  are  considered.  The  diversity  receiver  is  shown  in  Fig.4. 


t  =  T,2T,... 


r(t) 


►<&-►  \{.)dt   -^ 


i  + 


+ 


*■  (>y 


cos(2jrf  t) 


+1 


(.)' 


Hg>-* 


sin(2rf() 


max(Y2.„Y2.2 YZL| 

+ 

second-max  { Y2,i ,  Y2,2, . . . .  Y2>L  } 


choose 
largest 


max{Yu,Y1,2,....Y1.L} 
+ 
second-max{Yi_i,Y1-2,....YliL} 


{o, 


v, 


Figure  4.  Block  diagram  of  a  noncoherent  DPSK  receiver 

with  second  order  Post  -  Detection  Selection  Combining 
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A.        PROBABILITY  DENSITY  FUNCTIONS  OF  THE  DECISION  VARIABLES 

As  shown  in  Fig.  4  let  V,   and  V,   denote  the  decision  variables  for  the  signal 

branch  and  non-signal  branch,  respectively.  For  the  signal  branch,  let  V, ,  and  V,  2  denote 

the  two  dependent  random  variables  that  represent  the  two  largest  values  at  the  output  of 
the  L  noncoherent  DPSK  demodulators,  respectively 

Vhl  =max{Yl],Yl2....Y1L},  (66) 

Vj2  =  second  max  {^,,7,2....}^^).  (67) 

Also,  let  V2>1  and  V22  denote  the  two  dependent  random  variables  that  represent  the  two 
largest  values  at  the  output  of  the  non-signal  branch 

V2A  =  max{F21,F22....F2>J,  (68) 

V22  =  second  max {F21,F22....F2L  j.  (69) 

Then,  the  decision  variables  V,  and  V2  can  be  written  as 

V,  =  maxJFj , ,  Yu....YlL  }+  second  maxfc , ,  Yl2  ....Y1L  },  (70) 


or 


and 


or 


^1=^,1+^1,2,  (71> 


V2  =  max-jV^  j ,  Y22....Y2L  }+  second  max{F2 , ,  F2  2  ....F2X  ),  (72) 


Vi=Vu  +V2,2 ,  (73) 
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where  Ylk  and  Y2k ,  k  =  1,2,. ..L  are  the  amplitudes  at  the  outputs  of  the  DPSK  demodulator 
in  diversity  channel  k .  We  assume  that  Y]k  and  Y2k  are  independent  and  identically 
distributed  random  variables. 

For  the  signal  branch,  the  probability  density  function  of  the  k  th  signal  amplitude 

YXk  is  a  central  chi-square  density  function  given  in  (20) 


/«.(*»> -£T«P 


2^ 

2cr2 


***<>, 


(74) 


where  <72  is  given  in  (21)  and  (22)  as 


of  =4ol  +N0=2E  +  N0 


(75) 


V,  is  the  sum  of  Vu  and  V,  2 .  In  order  to  calculate  the  pdf  of  decision  variable  V, ,  we 


first  need  to  obtain  the  joint  pdf  of  V, ,  and  Vx  2  which  is  given  by  (Appendix  B) 
A.,  .vu  (vi-i '  v'.2  )  =  L(L-\)  fYu  (vu )  /yu  (vlt2 )  [Fyu  (v,  2  )f2 ,       v, ,  >  v1>2 , 


(76) 


where 


1,1 


1  (     v„  ^ 

At(vu)=— exp 


2cr,2 


VU^0' 


(77) 


/ru(vu)  = 


20"2 


exp 


'1.2 


2<j; 


v,-  >0, 


(78) 


and 


'*M-J 


Jo   2<Tf 


exp 


dyXk  =l-exp 


'     v      "i 
2?, 


(79) 
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The  cumulative  distribution  function  (cdf)  of  V,  is  found  by  integrating  the  joint  density 
fv   v  ,  (vi  i '  vi,2 )  over  tne  regi°n  shown  in  [3].  We  obtain 


h    "1-V1.2 


FV,M=   J        J/vu.Vu(VU»V1.2)  dvU  dvh 

L(L-l) 


0         v, 


'2    v,- 


=  1     J 


0         v, 


,  4^ 


•exp 


vu+vu 


2a? 


(     „..  \ 


-ii-2 


1-exp 


1,2 


2a; 


v        '  yj 


/<?>"  vi.i  ^  vi,2 


dvudvl2. 


(80) 


Performing  this  integration  yields  [3] 


FVi(v1)  =  L(L-l> 


1- 


A  -vv 


1+- 


2(7, 


2fff 


Z.-2 

+Z 


n-2^1 


(-1)^,0) 


where 


2  +  7      7 


r 


7(2+7) 


A   eXP 


V 


'.(2+7)' 

4a,2 


(81) 


(82) 


dFv(vx) 

Differentiating  (81)  with  respect  to  v, ,  fv  (v, )  = j—  - ,  yields  the  pdf  of  the  decision 


dvx 


variable  for  the  signal  branch  as  follows  [3] 


exp 


A(v,)=L(L-l)- 


t     2C7,2 


2cr,2 


L-l 


+ 


*rf     ? 


j  > 


(-iy 


1  — exp 


v,  7  * 


4*>    JJJ 


(83) 


Similarly,  for  the  non-signal  branch,  the  probability  density  function  of  the  k-th  noise 
amplitude  Y2k  is  a  central  chi-square  density  function  given  in  (25) 


fY"(yu)  =  2^&XP 


2a2 


y2k^o, 


(84) 
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where  a2  =Gn  =  N0 .  The  joint  pdf  of  V2>,  and  V22  is  given  by  (Appendix  B) 

/v2, ,v2.2  (vi,i .  v2,2 )  =  L (L - 1) /r2i  (v2 , )  /^  (v22  )  [fK2i  (v2>2 )}~2 ,      v2A  >  v22 ,         (85) 


where 


(  V  ^ 


44(v2.i)  =  — TexP 


2<t; 


1 


2,1 


2cr,2 


At(y2,2)  =  ^rexp 


K2.2 

2c722; 


v2,>0, 


v22>0, 


(86) 


(87) 


and 


v2.2  j 

FrJV2,2)=  J  T-Texp 


r  y 


)  2C2 


2* 


>\ 


2a227 


^     „       A 


4y2*  =l-exp 


2cr : 


v  z  J 


(88) 


To  find  the  pdf  of  V2  =  V2 ,  +  V2  2 ,  we  must  first  obtain  the  cdf  of  V2 .  By  integrating  the 
joint  density  fv   v  ^  (v2 , ,  v22 )  over  the  region  shown  in  [3]  we  obtain 


/2    v2-v22 


FV2  (V2  )  =    J         J  fvu  ,VM  (V2.1  .  V2.2  )  dvi.x  dv2.: 


for  v, .  >  v 


0  v2.2 


'2      V2-v2.2 


=  1       I 


L(L-l) 

,       4        eXP 

4(7, 


V2,l  +V2,2   ^ 


2cr2 


/L 


r   v    N\ 


-iL-2 


1-exp 


2.2 


V     2ffWJ 


2,1   ■"   K2.2 


dv2Adv2<2, 


(89) 


Performing  this  integration  yields  [3] 


F„(v2)=L(L-l) 


1- 


A       _V2/ 


1+- 


2cr 


V      —  2  y 


'2o 


L-^L-2^ 


/t=i 


+z  :  (-i)*w2(t) ,    (W) 


where 
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w2(k)=— — !«"%*  + 


2  +  k     k 


k(2+k) 


exp 


v2(2  +  k) 


4cr; 


(91) 


dFv(v2) 

Differentiating  (90)  with  respect  to    v2,    fv  (v2)  = ,  yields  the  pdf  of  the 

dvr, 


decision  variable  for  the  non-signal  branch  as  follows  [3] 


(      ,.     "\ 


exp 


fv(v2)=L(L-l)- 


2a: 


2a: 


&fL-2y-if( 


Aa22     f^[   k     )     k 


exp 


(92) 


B.       BIT  ERROR  PROBABILITY 

Using  the  results  derived  in  Part  A,  the  probability  of  error  for  second-order  Post- 
Detection  Selection  Combining  is 


'.-/ 


J7v2(v2)<fv2 


fv,  (vi )  dv\  ■ 


(93) 


Substituting  (83)  and  (92)  into  (93)  yields 


r     ..    \ 


P,=\\L2(L-V 


exp 


2a; 


0  v, 


2a[ 


Aa 


■  + 


&(L-2\ 


1  7=1 


v    J     J 


(-iy 


1  — exp 


/     ,.    > 


exp 


2a; 


x- 


2cr: 


v2       «(L-2>!(-l)*fi       J    va*^ 


4cj2 


+ 


2         *=1 


* 


1-exp 


4cr2 


V  -   J)l 


dv2  dvx 


(94) 
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The  result  is  found  in  terms  of  o2  and  o2  as 


K=L>(L-l)2\ 


OxG2 


+ 


{2{<j;-+<j*J   4(<t,2+<t22): 


L-2 


(L-2Y-\)k 
k  k 


4o: 


2(o,2+o22)2      (2  +  ^)(o,2(2  +  Jt)+2o22)2 


i=\  [  j      )    j 


2      2 

g:o; 


2o2o2 


2(o2+o22)2      (2o2+o2(2+;))2 


+ 


o, 


2(o2+o22)     2o2+o2(2  +  ;) 


^^fL-2YL-2^(-i)j+k 

+  Zj2j  — ■ — 

fc=i  ;=i  {    k     ){  j      )     jk 


X 


2o: 


o,  +o0 


2o2+o2(2  +  y) 

4o2  4o2 

(2  +  fc)(2o2+02(2  +  A:))     (2  +  £)(o2(2  +  £)+o2(2  +  ./)). 


(95) 


By  rewriting  the  parameters  ox  and  o2  as 


.2         /i^2 


and 


a^  =4oz+Aro=2£  +  Nc 


o-22=o-;=Ar0, 


(96) 


(97) 


the  bit  error  probability  is  found  in  terms  of  the  energy  per  diversity  channel  as 
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Pe=L2(L-lf 


IE 

1°. 

2E_ 


+  1 


+  2 


+ 


1 

(2E     ^ 
—  +  2 


+lfL~2^ 


k=l 


\  J 


1 


f2E       V 

—  +  2 


(2  +  k) 


LV 


2£ 

Nn 


+  1 


(2  +  ik)+2 


,ff^-^(-iy 

J    J   ; 


7=1 


V     J        J 


2£ 


+  1 


<2E 


Nr 


+  1 


2£ 


+2  — +4+j 


+ 


1 


1 


v 


IE 
—  +  2 
N, 


—  +  4  +  j 
Nn 


j 


JJ 


^k^(L-2)(L-2)(-\)J+k 


k=\  j=\ 


k 
1 


v    J     J 


jk 


-  +  2     -  +  4  +  J     {2  +  k) 


+ 


2  + 


(2E       ^ 


yN0 


+  1 


{2  +  k) 


(2  +  k) 


(2E 


LV 


Nr 


+  1 


{2  +  k)+{2  +  j) 


(98) 


—        E* 
Substituting  E  by  in  (98),  the  bit  error  probability  becomes 

1-/ 
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pb=L2(L-iy 


( 

2 


2Eb_ 
LN0 

2Eb 


+  1 


+  2 


+  ■ 


2Eb 


+  2 


+  l(L~2] 


k=\ 


2Eh 
LNn 


+  2 


(2  +  k) 


L\ 


2Eb 
LNr 


+  1 


(2  +  k)+2 


+  Y2(L-2)(r}L 
m[    J     J 


J 


2Eh 

LNr 


+  1 


AEb 

LNr 


+  4  +  7 


AEb 

LNn 


+  4  +  y 


i=i  >=i  ^   *     X    J     )     jk 


1 


x 


2Eb      ^      AEb      A      .      ,        s 
+  2  -  +  4+;      (2  +  jfc) 


Ltfr 


+ 


ZJV, 


2  + 


2£fc 
ZJV\ 


+  1 


(2  +  *) 


(2  +  *:) 


2Eh 

LNr 


+  1 


(2  +  *)+ (2  + j) 


(99) 


In  the  next  chapter  we  will  derive  the  bit  error  probability  when  third  order  Post- 
Detection  Selection  Combining  is  employed. 
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V.       THIRD  ORDER  POST-DETECTION  SELECTION 
COMBINING  (PDSC-3)  ANALYSIS 

In  the  third-order  Post-Detection  Selection  Combining  (PDSC-3)  technique,  the 

three  largest  values  of  L  diversity  samples  at  the  outputs  of  the  L  noncoherent  DPSK 
demodulators  are  considered.  The  diversity  receiver  is  shown  in  Fig.5. 
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u- 
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+ 
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+ 

third-max{  Y2 1  ,Y2  2, . . .  Y2.i.} 
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+  1*         ma 
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Figure  5.  Block  diagram  of  a  noncoherent  DPSK  receiver 

with  third  order  Post  -  Detection  Selection  Combining 
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A.        PROBABILITY  DENSITY  FUNCTIONS  OF  THE  DECISION  VARIABLES 

As  shown  in  Fig.  5  let  V,  and  V2  denote  the  decision  variables  for  the  signal  branch  and 
non-signal  branch,  respectively,.  For  the  signal  branch,  let  Vtl,V12  andV13  denote   the  three 

dependent  random  variables  that  represent  the  three  largest  values  at  the  output  of  the  L 
noncoherent  DPSK  demodulators,  respectively 

Vu  =max{ril,712....ri  J,  (100) 

Vl2  =  second  max{F11,ri2....F1  L\  (101) 

Vu  =  third  max{ru,yi2....yu}.  (102) 

Also,  let  V2l ,  V2 ,  and  V23  denote  the  three  dependent  random  variables  that  represent  the  three 

largest  values  at  the  output  of  the  non-signal  branch 

V2A  =max{y2il,F2i2....r2,J,  (103) 

V22  =  second  max{y2],F22....F2X  j,  (104) 

V23  =  third  max{r21,y22....r2t}.  (105) 
Then,  the  decision  variables  V,  and  V2  can  be  written  as 

V,  =max{yu,yi2....yiL}+ second  max{yu,FI2....F1L}+  third  maxj^^^^....^^},  (106) 


or 


v;=vu+y1>2+y1(3,  (107) 

and 

V2  =max{l/21,F22....F2jL}+ second  max{F2],y22....y2L}+  third  max{y2],y22....y2X},      (108) 


or 


v2  =V1A  +V2,2  +V23,  (109) 
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where  YXk  and  Y2k ,  k  =  1,2,. ..L  are  the  amplitudes  at  the  outputs  of  the  DPSK  demodulator  in 
diversity  channel  k .  We  assume  that  Ylk  's  and  Y2k '  s  are  independent  and  identically 
distributed  random  variables. 


For  the  signal  branch,  the  probability  density  function  of  the  k-th  signal  amplitude  F, 


u 


is  a  central  chi-square  density  function  given  in  (20) 


2o\ 


2d,2 


y**o, 


(110) 


where  <72  is  given  in  (21)  and  (22)  as 


<72=4cr2+Aro=2£  +  AV 


(111) 


Since  V]  =  V, ,  +  Vl2  +  Vu ,  we  need  to  obtain  the  joint  pdf  of   Vu ,  VX2  and  Vj>3  which  is  given 
by  (Appendix  C) 


/vu,vu,vI.3(vu'vi,2.v1>3)  =  L(L-l)(L-2)/yu  (vu)/yu  (v1>2)/Fu  (v13)  [fYu  (vu)}~3 , 


Vl.l   ^  V1.2  ^  V1.3  . 


(112) 


where 


/ru(vu)=  — rexp 

2(7," 


f 


1.1 


2(7, 


vu'*0, 


(113) 


■Mvu)  =  — TexP 


1,2 


2cr2 


v12>0, 


(114) 


/y"(Vu)=2^TeXP 


1.3 


2cr; 


v,,>0, 


(115) 


and 


i 


1_ 


2C72 


rfyu  =  1-exp 


1,3 


2a; 


(116) 
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To  find  the  pdf  of  Vx ,  we  must  first  obtain  the  cumulative  distribution  function  (cdf)  of  V, .  By 
integrating  the  joint  density  fv   v   v   (vltl ,  vl2 ,  v13 )  over  the  region  where  Vj ,  >  Vj  2  >  v,  3  and 

Vj  j  +  v,  2  +  v,  3  <  v,  as  shown  in  [3]  we  obtain 


v\/  (vi-vu)/ 

/3  /2  vi-vi.2-vi.3 

FV,(Vl)=J  J  J      A,,V,2.VU(V1,1'V1.2.V1,3)^U^1,2^1.3 

0  "1.3  v1.2 


(117) 


Performing     the     integration     and     differentiating     the     result     with     respect     to     v, , 

dFv{vx) 

fv  (vi ) '  the  pdf  of  the  decision  variable  for  the  signal  branch  is  given  by  [3] 

1  dv, 


,    ,L{L-\\L-2)e2^ 

•Mlj""  2  "2d,2 


kZ(L-3\ 


W*\  )2 


+ 


Jt=l 


(-1)' 


v        ; 


fc2 


fcv, 
2cr2 


-3 


1-exp 


6d2, 


(118) 


Similarly,  for  the  non-signal  branch,  the  probability  density  function  of  the  k-th  noise 
amplitude  Y2k  is  a  central  chi-square  density  function  given  in  (25) 


/r2i(3'2jfc)  =  T-TexP 

2<T, 


r   v    ^ 

2C72 


y2^o, 


(119) 


i         ~i 


where  cr^  =  G~n  =  N0  .  The  joint  pdf  of  V2 ,  ,  V,  2  and  V23  is  given  by  (Appendix  C) 
fv2Av2a,v2^v2A>v2,2'v2,3)  =  L(L-l)(L-2)fYu(v2tl)fy^^ 

V2A   ^  V22  >  V23, 

where 


(120) 
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/y2,(v2,,)  =  T-^eXP 


fy2t  (v2.2) 


2a: 


2a: 


2c\t 


exp 


(    v      > 
"2,2 


2(j; 


v    ~ 2  ) 


fr    (v2s)  = 7exP 

u  2a\ 


(     v      ^ 
^2,3 


2cr; 


v2,  >0, 


v22>0, 


v2,3>0, 


(121) 


(122) 


(123) 


v         z  J 


and 


1 


1_ 

)     2fT2 


dy2k  =l-exp 


(     v      ^ 

1/2,3 


2a22; 


(124) 


To  find  the  pdf  of  V2 ,  we  must  first  obtain  the  cdf  of  V2 .  Integration  of  the  joint  density 
fv2Uv21>v2Av2,i ^2,2^2,3)  over  the  region  where  v21  >v22  >v23  and  v21  +v22  +v23  <v2  [3] 


leads  to 


v2/     W2  -v2 


,,_)/ 


V2-V2f2-V2t3 


FV2(V2)  =    J  J  J      /v2,,V2.2.V2.3(V2,l'V2,2'V2,3)^V2,l^V2,2^V2,3 

0  v2.3  v2.2 


(125) 


Performing     the     integration     and     differentiating     the     result     with     respect     to      v2, 

d  Fv  (v2) 

/v  (v2 )  = = — — ,  the  pdf  of  the  decision  variable  for  the  non-signal  branch  is  given  by  [3] 

2  dv-, 


fvM)= 


L{L-l\L-2)e  2g| 

2  "  2cr22 


x^ 


^YL-3V-iy 

v  7    y 


6(2C722)2        ;=l 


;' 


2a: 


-3 


1-exp 


jv2 


(126) 
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B.       BIT  ERROR  PROBABILITY 

Using  the  results  derived  in  Part  A,  the  probability  of  error  for  third-order  Post- 
Detection  Selection  Combining  is 


i-J 


j7v2(v2)^ 


/v,  (vi )  dv\ 


(127) 


Substituting  (118)  and  (126)  into  (127)  yields 


Pe=}\L\L-Xf{L-2f 


0  v, 


-    2    ,exp 

16<7,(72 


2cr2     2(72  ; 


x< 


*■< 


6(2e?J 


f      +Wi-3U 


*=1 


£-3, 


?+I 


L-3 


6(2(j22)2     £(    J 


(-iy 


2cr,2 


JV2 

2(72 


-3 


-3 


l-e 


l-e 


6cr,2 


/J 


jv2     "\ 


6ct4 


/J 


dv-,dv 


VA,V\ 


(128) 


The  probability  of  error  is  found  in  terms  of  cfanda2  as 


P,  =L2(L-l)2(L-2f[A  +  B  +  C  +  D], 


(129) 


where 


A  = 


_4_6 
(7,  (T2 


,-r2-r6 

(7,072 


6(c72+a2)5  +12(c72+C72)4  +  36{a2 +<J22J  ' 


Z.-3 


(L-3M-1) 


k  =  \ 


k     J   £2 


fc<74<74 


4       2 
<7,  (7, 


27(7,4C72 


4(<7,2  +  a2 )4     4(<7 2  +  a2  /      4(3(7 2  +  3(72  +  ko\  f 


+ 


ko]o\ 


(7,(72 


9C72C72 


6((7 2  +  cr22  /      4((72  +  crj  J       4(3(7 ,2  +  3(7  2  +  &<72  j 


3(7 : 


K(72  (72  3(72 

+ 12((72  +(72 )2  " 4(c72  la; )+  4(3(72  +  3(72  +£(72 ) 
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c.w«™ 


J 


J 


i°l 


4(a2  +  a2)4  +12(a2  +<j2J     4((T,2  +a2)3 


+ 


81a, 


4(3  +  7'X3^2+3a22+ya12)3 


and 


D^^(L-3YL-3"\(-lf+j 


k=l  j=\ 


A 


J 


w 


kjG\c\ 


3jG2a2 


+  ■ 


27  jc2a22 


2(a2  +  a2  J     4(a2  +  a2  J     4(3a,2  +  3a22  +  ka\  J 


+ 


kjo\  3ja22       ,  9y'a2 


4(a2+a22)2     4(a2  +  a2)    4(3a2  +  3a2  +  J:a22) 


3£a' 


+ 


9a: 


27a: 


4(a2+a2)2     4(af+a2)    4(3a2  +  3a2+£a2) 

81fca2 81a2 

+  4(3  +  7X3a2  +  3a2  +  ;a2J    4(3  +  y)(3a2  +  3a2  +  ja2) 

81a2 

4(3  +  j  X3a2  +  3a2  +  jo2  +  ka\ ) 


+ 


Rewriting  the  parameters  a,  and  a2  as 


o-2=4a2+7V0=2£  +  Ar0=2-^  +  Ar0  , 


(130) 


and 


c72=a2=7V0, 


(131) 


allows  an  expression  for  the  bit  error  probability  in  terms  of  Ebl  N0  as 
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pb  = 


L2(L-l)2(L-2)2  f   1    3a2+3ab  +  2b2     &(L-3\{-\)k 
4  ]  144  b5  +^l     k 


k=\ 


X 


1   3ka(a  +  \)+2b(2ka-3a2+k-3)+4b2(k-3a-6b)  81 
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+ 


+ 


L-3  L-3 

II 

*=1   ;=1 


L-3^ 


'1-3^ 


v    ■/    ; 


(-D 


3(9a 2  +  3a  (6Z7  +  k)+  (6b +kf ) 

(6/?  +  kf 
'3aj  +  (6b  +  k)(j-3)     kja  +  (j-3)kb-3jab  +  (9-3j)2b: 

y : — r- 


(3  +  k)(a(3  +  k)+3)2 


t+j 


w 


4b : 


+  81 


(6b  +  kf 

k-((3+j)a  +  3)      | - 

(3  +  j)((3  +  j)a  +  3f      (3  +  jfa  +  (3  +  j)(3  +  k)j[ 


(132) 


where  a  = 


(2I      ^ 


ZJV„ 


+  1 


and£  = 


In  the  next  chapter  we  will  present  the  performance  of  the  various  diversity  combining 
techniques  previously  discussed. 
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VI.      NUMERICAL  RESULTS 

All  expressions  related  to  the  bit  error  probability  and  probability  density 
functions  of  the  decision  variables  are  analytically  derived.  Numerical  results  are 
obtained  by  using  MATLAB  [13].  A  signal-to-noise  ratio  range  of  6-20  dB  for  is  used  in 
the  numerical  evaluations. 

A.        COMBINING  TECHNIQUES  AND  DIVERSITY 

We  obtained  numerical  results  for  Post-Detection  Selection  Combining,  Selection 
Combining  and  Equal  Gain  Combining.  As  discussed  in  the  preceding  chapters,  the  main 
subjects  of  interest  are  :  first,  second  and  third  order  Post-Detection  Selection  Combining. 

The  test  cases  allow  us  to  compare  the  combining  techniques  for  different  orders 
of  diversity,  L.  Plots  for  Post-Detection  Selection  Combining,  Selection  Combining  and 
Equal  Gain  Combining  for  different  orders  of  diversity  are  given  in  Fig.6  through  Fig.  12. 
The  bit  error  probability  expressions  given  in  [3]  are  used  for  the  numerical  evaluation  of 
Selection  Combining. 

To  see  the  differences  in  the  performance  among  these  three  techniques,  we 
compare  PDSC-1,  SC-1,  and  EGC  for  diversity  order  L=l,  2,  3,  4  and  5  in  Fig.  13  through 
Fig.  17.  Then  we  compare  PDSC-2,  SC-2  and  EGC  performance  for  L=2,  3,  4,  5  and  6  in 
Fig.  18  through  Fig. 22. 

Finally,  PDSC-3,  SC-3  and  EGC  performance  are  compared  for  L=3,  4,  5,  6,  7,  8, 
9  and  10  respectively  in  Fig.23  through  Fig. 30.  In  addition  to  L=3,  4,  5,  6,  7,  diversity 
orders  of  8,  9  and  10  are  employed  to  show  a  slight  improvement. 
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B.      COMPARISONS 

We  note  that  EGC  performs  similarly  or  better  than  the  other  two  methods. 
However,  our  main  objective  is  to  show  the  performance  difference  between  PDSC  and 
SC. 

In  Fig.6  through  Fig.  12,  performance  curves  of  PDSC- 1,  PDSC-2,  PDSC-3,  SC-1, 
SC-2,  SC-3  and  EGC  are  plotted  for  different  values  of  diversity  orders.  It  is  seen  that  as 
the  diversity  order  increases,  the  performance  of  the  receiver  increases. 

In  Fig.  13  through  Fig.  17,  PDSC-1  and  SC-1  allow  comparison  for  L=l,  2,  3,  4  and 
5.  When  L=l,  (i.e.  when  there  is  no  diversity),  all  three  methods  perform  the  same  as 
seen  in  Fig.  13.  Starting  from  diversity  of  L=2,  we  notice  a  slight  performance  difference 
between  PDSC-1  and  SC-1,  PDSC-1  performs  better  than  SC-1,  and  worse  than  EGC  as 
previously  stated.  As  the  value  of  L  increases  the  performance  difference  between  PDSC- 
1  and  SC-1  increases  in  favor  of  PDSC-1  as  shown  in  Fig. 14,  15,  16  and  17. 

The  performance  of  PDSC-2,  SC-2  and  EGC  is  illustrated  in  Fig.  18  through 
Fig. 22.  Since  two  signals  are  combined,  these  figures  are  plotted  for  diversities  of  2 
through  6. 

When  L=2,  as  seen  in  Fig.  18,  all  combining  techniques  perform  the  same  as 
expected.  For  L=3,  there  is  virtually  no  difference  between  the  three  techniques.  In  Fig. 
20  (i.e.,  L=4),  and  in  Fig.21  (i.e.,  L=5),  we  see  that  PDSC-2  performs  essentially  the 
same  as  SC-2.  Again,  EGC  performs  better  than  the  other  two  methods.  In  Fig.22 
(i.e.,L=6),  PDSC-2  is  slightly  better  than  SC-2. 

The  performance  for  SC-3,  EGC  and  PDSC-3  is  shown  in  Fig. 23  through  Fig.30. 
In  order  to  more  clearly  see  the  performance  differences,  plots  are  given  for  diversities  of 
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L=3,  4,  5,  6,  7,  8,  9  and  10.  Since  three  signals  are  combined  for  PDSC-3  and  SC-3,  the 
performance  curves  are  plotted  starting  with  L=3.  For  L=3,  all  three  methods  perform 
exactly  the  same.  Although,  it  is  really  difficult  to  distinguish  the  performance  curves  in 
Fig.24-26,  very  slight  differences  in  the  curves  indicate  that  EGC  performs  better  than 
SC-3,  and  SC-3  performs  slightly  better  than  PDSC-3. 

Beginning  with  Fig.27,  evaluated  for  L=7,  we  see  that  PDSC-3  is  starting  to 
perform  slightly  better  than  SC-3  for  an  SNR  value  of  17  dB.  As  the  diversity  orders  are 
increased  to  L=8,  9  and  10  the  performance  differences  between  PDSC-3  and  SC-3 
become  more  significant,  favoring  PDSC-3.  For  L=10,  PDSC-3  outperforms  SC-3  for 
SNR  values  of  12  dB  or  larger.  But  the  difference  is  small. 

In  summary,  PDSC-1  performs  better  than  SC-1,  while  PDSC-2  and  PDSC-3 
perform  essentially  the  same  as  SC-2  and  SC-3,  respectively.  EGC,  as  expected,  performs 
better  than  the  other  two  techniques.  For  time  diversity,  the  EGC,  PDSC  and  SC  methods 
require  only  one  receiver.  On  the  other  hand,  for  space  and  frequency  diversities  EGC 
and  PDSC  methods  require  L  receivers.  Although  the  SC  method  requires  only  one 
receiver  for  time,  space  and  frequency  diversities,  it  does  require  a  more  complex 
receiver,  one  which  has  the  pre-detection  combiner. 
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Figure  6.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading  channel 
using  first  order  post-detection  selection  combining  (PDSC-1)  for  diversity  orders  of 
L=l,2,3,4and5. 
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Figure  7.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading  channel 
using  second  order  post-detection  selection  combining  (PDSC-2)  for  diversity  orders  of 
L=2,  3,  4,  5  and  6. 
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Figure  8.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading  channel 
using  third  order  post-detection  selection  combining  (PDSC-2)  for  diversity  orders  of 
L=3,  4,  5,  6  and  7. 
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Figure  9.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading  channel 
using  first  order  selection  combining  (SC-1)  for  diversity  orders  of  L=l,  2,  3,  4  and  5. 
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Figure  10.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading 
channel  using  second  order  selection  combining(SC-2)  for  diversity  orders  of  L=2,  3,  4,  5 
and  6. 
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Figure  11.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading 
channel  using  third  order  selection  combining  (SC-3)  for  diversity  orders  of  L=3,  4,  5,  6 
and  7. 
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Figure  12.  Performance  of  the  noncoherent  DPSK  receiver  over  a  Rayleigh  fading 
channel  using  equal  gain  combining  (EGC)  for  diversity  orders  of  L=l,  2,  3,  4  and  5. 
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Figure  13.  Receiver  performance  of  PDSC-1,  SC-1,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=l. 
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Figure  14.  Receiver  performance  of  PDSC-1,  SC-1,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=2. 
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Figure  15.  Receiver  performance  of  PDSC-1,  SC-1,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=3. 
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Figure  16.  Receiver  performance  of  PDSC-1,  SC-1,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=4. 
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Figure  17.  Receiver  performance  of  PDSC-1,  SC-1,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=5. 
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Figure  18.  Receiver  performance  of  PDSC-2,  SC-2,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=2. 
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Figure  19.  Receiver  performance  of  PDSC-2,  SC-2,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=3. 
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Figure  20.  Receiver  performance  of  PDSC-2,  SC-2,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=4. 
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Figure  21.  Receiver  performance  of  PDSC-2,  SC-2,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=5. 
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Figure  22.  Receiver  performance  of  PDSC-2,  SC-2,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=6. 
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Figure  23.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=3. 
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Figure  24.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=4. 
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Figure  25.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=5. 
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Figure  26.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=6. 
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Figure  27.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=7. 
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Figure  28.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=8. 
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Figure  29.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=9. 
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Figure  30.  Receiver  performance  of  PDSC-3,  SC-3,  and  EGC  over  a  Rayleigh 
fading  channel  for  L=10. 
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VII.    CONCLUSIONS 

The  main  contribution  of  this  thesis  is  the  analysis  of  the  new  Post-Detection 
Selection  Combining  Technique.  The  expressions  for  the  bit  error  probability  are  derived 
for  First,  Second  and  Third  order  Post-Detection  Selection  Combining  DPSK  receivers. 
Numerical  Results  are  compared  with  the  results  obtained  from  Selection  Combining  and 
Equal  Gain  Combining.  It  is  observed  that  combining  the  two  (or  three)  largest  signals 
(PDSC-2,  PDSC-3,  SC-2  or  SC-3  )  offers  significant  performance  improvement  over 
selecting  the  maximum  signal  (i.e.  PDSC-1,  SC-1). 

First  order  Post-Detection  Selection  Combining  (PDSC-1)  which  is  done  by 
selecting  the  largest  signal  performs  better  than  First  order  Selection  Combining  (SC-1). 
As  the  order  of  the  diversity  increases,  PDSC-1  outperforms  SC-1.  The  improvement 
increases  as  the  diversity  increases. 

Second  order  diversity  combining  is  performed  by  combining  the  two  largest 
signals  as  PDSC-2  and  SC-2.  Numerical  results  indicate  that  SC-2  performs  the  same  as 
PDSC-2  for  diversity  values  of  L=2  to  4.  However  beginning  L=5,  PDSC-2  offers 
improvement  over  the  performance  of  SC-2.  We  can  conclude  that  as  the  order  of 
diversity  increases  performance  improvement  of  PDSC-2  becomes  more  significant 
relative  to  the  SC-2.  Again,  PDSC-2  performs  only  slightly  worse  than  the  Equal  Gain 
Combining. 

In  the  PDSC-3  and  SC-3  case,  numerical  results  show  that  the  SC-3  performs 
essentially  the  same  as  PDSC-3  for  all  diversity  orders.  Both  PDSC-3  and  SC-3 
performances  are  slightly  worse  than  EGC. 


69 


In  summary,  Post-Detection  Selection  Combining  and  Selection  Combining 
receivers  perform  worse  than  EGC  receivers  (up  to  1.5  dB  in  signal-to-noise  ratio).  On 
the  other  hand,  selection  combining  techniques  are  independent  of  the  diversity  order  L. 
This  is  desirable  because  L  may  vary  with  location,  as  well  as  time.  For  example,  in 
practice  different  diversity  channels  may  have  different  path  lengths.  Thus  the  signal-to- 
noise  ratio  of  each  diversity  channel  can  be  different  resulting  in  unequal  gain  combining. 
This  can  cause  a  serious  degradation  for  EGC  receivers.  SC-1  and  PDSC-1,  SC-2  and 
PDSC-2,  and  SC-3  and  PDSC-3  have  a  robust  performance  as  long  as  L  is  greater  or 
equal  to  1,  2  and  3,  respectively.  However,  the  SC  method  requires  a  more  complex 
receiver  than  the  EGC  and  PDSC  methods,  since  it  requires  a  pre-detction  combiner. 
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APPENDIX  A 

DERIVATION  OF  THE  PROBABILITY  DENSITY  FUNCTION  FOR  THE 
LARGEST  RANDOM  VARIABLE 

In  this  appendix  ,  the  pdf  expressions  are  derived  for  largest  variable  in  a  set  of  L 
i.i.d.  random  variables.  We  assume  that  Xl,X2,X3,...,XL  are  L  independent,  identically 
distributed  random  variables  with  probability  density  functions 

fX]  (*i  )=  fx2  fe)=  fx3  fa )=  -  =  fxL  (xl)=/x  W-  (A.l) 

An  expression  will  be  derived  for  the  largest  random  variable  in  a  set  of  L  i.i.d. 
random  variables  as  given  by 

Z  =  max{X„X2,X3,...Xj.  (A.2) 

In  terms  of  order  statistics,  it  can  be  interpreted  as  [1 1] 

Xrig)<Zr2g)<Xrjg)<...<Xrtg)<...<Xrig).  (A3) 

We  next  form  the  L  random  variables  z,  such  that 

Zi  =  *r,  ^z2  =  xr2  < <  z,  =  xTt  < <zL=  xrL  (A.4) 

The  k-th  largest  random  variable  in  this  sequence  is  called  the  k-th  order  statistic.  We 
determine  its  density  function  by  using 

fz{z)dz  =  ?r{z<zk<z  +  dz}.  (A.5) 

This  expression  simply  describes  an  event  that  occurs  when  k-1  of  the  random  variables 
are  less  than  z,  and  the  L-k  random  variables  are  larger  than  z+dz.  In  other  words,  we 
have  three  subsets  with  associated  probabilities  as  follows: 
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D}  =  {x  <  z} >  Px  =  Prfo  )=Fx(z) 

D2={z<x<z  +  dz} >  p2  =  Pr(Z)2 )  =  fx  {z)dz  (A.6) 

D3  ={z  +  dz<x} >  p,  =Pr(D3)  =  l-  Fx{z  +  Az), 

where  Fx  (jc)  is  the  cumulative  distribution  and  /x  (jc)  is  the  density  function  of  x.  We 
can  define  the  order  statistics  "Dl  occurs  fc.  times"  as  [1 1] 


P\  '  Pi  2  P3 3     '    ^1  +  ^2  +  ^3  =  ^  •  (A.7) 


kx\k2%\ 


For  the  event  defined  above,  D,  occurs  k-1  times,  D2  occurs  once,  and  Z)3  occurs  L-k 

times,  hence 

kx=k-\       k2=\       k3=L-k.  (A.8) 

Substituting  these  expressions  in  (A.7),  we  obtain 


fz(z)dy  =  ?r{z<zk<z  +  dz}=  tk_xyl_k  yFxk''(z)fx{z)dz[l-Fx(z  +  dz)]L-k.(A.9) 


As  a  general  expression  ,  the  pdf  of  the  k-th  order  statistic  is  obtained  as 

In  this  case,  the  largest  value  corresponds  to  k=  L  and  the  probability  density 
function  of  the  largest  random  variable  in  L  random  variables  is  found  to  be 


fz(z)=Lfx{z)FxL-\z).  (A.11) 
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APPENDIX  B 

DERIVATION  OF  THE  JOINT  PROBABILITY  DENSITY  FUNCTION  FOR 
THE  LARGEST  AND  THE  SECOND  LARGEST  RANDOM  VARIABLES 


This  time  we  are  looking  for  the  joint  probability  density  function  of  the  largest  and 
the  second  largest  random  variable  in  a  set  of  L  random  variables 

Z,  =  first  max{X, ,  X2 ,  X3  ,...XL}, 
Z2  =  second  max{Xl,X2,X3,...X  L}. 
The  probability  density  function  of  Zj  is  known  from  the  previous  part  as 

/z.feJ^/xfe.K^fo)-  (B.l) 

For  Z2  ,  k  =  L  - 1  and  the  probability  density  function  of  Z2  is  obtained  by  substituting 
k  =  L-\  in(A.7) 


fz  fe)=77 x     w      7 rrFxL-2(z2)[l-Fx(z2)]L-{L-1]  fx(z2).  (B.2) 

2  ((L-l)-l)!(L-(L-l))! 


Simplifying  the  expression  above,  we  obtain 

fzM)=L(L-\)fx{z2)FtL-2{z2)[\-Fx{z2)].  (B.3) 

To  find  the  joint  pdf  for  Zi  and  Z2  ,we  need  to  define  a  series  of  events  as  follows 

fztj^(zi*z2)dz1dz2  =Pr{z!  <Z,  ^z.+Jz^z,  <Z2  <z2+*2},  (B.4) 
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Dx  =  {x<z2} >/>,  =Fx(z2) 

D2  ={z2  <x<z2+dz2} >p2  =/,(z2)<*Z2 

D3={z2+dz2<x<z1} >P3=Fx(zi)-Fx(z2+dz2)  (B.5) 

D4={Zl<x<  z,  +dz, } >  p4  =  /,  (z,  )dzt 

D5  ={x>zl+dzl} >p5  =l-Fx(zl+dz1). 


In  this  case,  Dj  occurs  (L-2)  times  and  ki=  L-2  ,  D2  occurs  once  and  k2=l,  D3  does  not 
occur  and  k3=0,  D4  occurs  once  and  k4=l,  D5  does  not  occur  and  ks=0.  Using 


Prfc,*a>...0=        L!        p*  p2k2  ....p3kr     ,     kx+k2....  +  kr  =L,  (B.6) 

we  obtain 


/z1,22fel^2K1*2=Prfel   <Z1   <Z,+*1,Z2<Z2   <Z2+JZ2} 

=  (L_2)!  F^2  (Z2  )A  (Z2  )*2  A  fc  )*1  (B  ?) 

=  l(L- l)/x (z2 )/,  (z,  )F^2 (z2  )rfz1dfe2 . 


Thus 


fzuz2  fc ,  z2 )=  L(L-l)/,  (z,  )/x  (z2  )F^2  (z2 )  •  (B.8) 
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APPENDIX  C 

DERIVATION  OF  THE  JOINT  PROBABILITY  DENSITY  FUNCTION  FOR 
THE  FIRST,  SECOND  AND  THIRD  LARGEST  RANDOM  VARIABLES 

In  the  final  case,  we  want  to  find  an  expression  for  the  joint  pdf  of  the  first,  second 
and  third  largest  random  variables  in  a  set  of  L  random  variables. 

Zj  =  first  max{Xl,X2,X3,...X  L}, 
Z2  =  second  ma.x{Xl,X2,X3,...X L}, 
Z3  =  third  max{Xj  ,X2,  X3,...XL  }. 

Similar  to  the  other  cases,  the  joint  pdf  of  Z\ ,  Z2  and  Z3  can  be  defined  as 

/z1,z2,z3fc^2,z3)=Pr{z,  <Z,  <zx+dz,,z2  <Z2  <z2+dz2,z3  <Z3  <z3+dz3}, 

and  in  terms  of  sets  of  events  as 

D1={x<z3} >P1=Fx(z3) 

D2={z3<x<z3  +  dz3 } >  p2  =  fx  (z3 )  dz3 

D3={z,+dz3<x<z2} >Pi=Fx(z2)-Fx(z3+dz3) 

D4={z2<x<z2+dz2} >P,=fx{zi)dz2  (C.l) 

D5={z2+dz2<x<Zl} >p5=Fx(z1)-Fx{z2+dz2) 

D,={zx{x<zx+dzx} >P6=fx(zi)dzl 

D1={x>z1+dzl} >p7=l-Fx(zl+dz1)    . 

In  this  particular  case,  Di  occurs  (L-3)  times  and  ki=  L-3,  D2  occurs  once  and 
k2=l,  D3  does  not  occur  and  k3=0,  D4  occurs  once  and  k4=l,  D5  does  not  occur  and 
k5=0,  D6  occurs  once  and  k6=l,  D7  does  not  occur  and  k7=0  . 
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Using  (B.6)  we  obtain 
fzl,z„Azl,z2,z3)  =  L(L-lXL-2)fx(z1)fx{z2)fx(z^FxL-3(z3).  (C.2) 
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